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Abstract 



In the core of the vortex of a superconductor, energy levels appear inside the 
gap. We discuss here through a random matrix approach how these levels are 
broadened by impurities. It is first shown that the level statistics is governed 
q ' by an ensemble consisting of a symplectic random potential added to a non- 

random matrix. A generalization of previous work on the unitary ensemble 
in the presence of an external source (which relied on the Itzykson-Zuber in- 

r , tegral) is discussed for this symplectic case through the formalism introduced 

C*~) | by Harish-Chandra and Duistermaat-Heckman. This leads to explicit formu- 

~j lae for the density of states and for the correlation functions, which describe 

(^ • the cross-over from the clean to the dirty limits. 
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I. INTRODUCTION 

The energy levels inside a vortex of a superconductor have been characterized long ago 
KJ, but recent studies have dealt with the broadening of these levels by impurities. The 
density of states has been computed both in the clean limit and for the dirty case. 

In this article, we discuss the crossover between these two limits. The impurity potential 
is handled through a random matrix theory. However the matrix elements between different 
energy levels due to the impurities are strongly correlated, and therefore one is far from a 
usual Wigner-Dyson random matrix theory. 

It is shown below that the crossover is characterized in this case by a random symplectic 
matrix coupled to an external non-random source matrix. For handling this problem we 
develop a technique which generalizes earlier work on the unitary ensemble in which one 
considered the cross-over from a deterministic hamiltonian to a fully random hermitian 
hamiltonian 0-§[. In the hermitian case, the formalism relied on an integral over the 
unitary group due to Harish- Chandra || and rederived in the context of random matrix 
theory by Itzykson-Zuber J7|. All the correlation functions for the energy levels were then 
found explicitely. For the vortex problem considered here the perturbation due to impurities 
is a matrix with a symplectic structure, which we treat as a random potential added to the 
unperturbed diagonal matrix consisting of the regularly spaced energy levels. It is shown 
below that for this problem a similar integral over the symplectic group, considered by 
Harish- Chandra ||, and more recently generalized by Duistermaat and Heckman || leads 
also to explicit formulae. The crossover from the clean spectrum to the dirty limit follows 
from this formalism. We will show how this general behavior for the crossover is relevant to 
the actual problem of the excitation energy spectrum in a superconducting vortex. 

II. SYMPLECTIC STRUCTURE OF THE PERTURBATION FOR THE ENERGY 

LEVELS INSIDE A VORTEX 

The energy eigenvalues and eigenfunctions of a quasi-particle inside the vortex of a su- 
perconductor, were obtained long ago |IJ . The two component wave functions of excitations 
■0 = (01, fa) in a superconductor satifies the Bogolubov-de Gennes equation, 

[a z {H + Vimp) + a x ReA{r) + a y ImA(r)]i> = E^ (2.1) 

where H = p 2 /2m — Ep, and A(r) is a gap order parameter. The impurity potential Vi mp 
is a sum of short range scattering sites rf. 

V imp (r) = ^2V t 5(r-r l ) (2.2) 

i 

In the absence of impurity, the Schrodinger equation ( [2.1|) leads to a spectrum of equidistant 
states E® 

E° n = -hu (n - i) (2.3) 



where n is an integer defining the quantized angular momentum of the vortex, and Kujq = 
A 2 /Ep. We assume that the spacing hu between these levels is much smaller than the gap 
A, and consequently there are many excitation levels in the vortex core. 

Using the explicit wave functions for the unperturbed eigenstates, one finds that the 
matrix elements of the interaction due to the i-th impurity are 

y -2K( n ) 

Km = - ; l e^ m ~ n ^[Jn(k F r t )J m (kFr^ - J n -i(k F ri) J m _i(A; F r,)] (2.4) 

Apt; 

where (r*, (pi) is the position of the i-th impurity in polar coordinates. The indices n and m 
run from — N to +iV where iV ~ A/uoq ~ Ep/ A is assumed to be a large number. 

The scattering time r, given by 1/r = 2niV 2 m, where rii = iVj/£ 2 is the density of 
impurities in the vortex ; JVj the total of impurities in the cross-section of a vortex of radius 

Three different situations may occur: a) a dirty limit A << 1/r << Ep, b) a clean limit 
ujq < 1/r << A, c) a superclean case 1/r << u Q ~ A 2 /E F . However, in ]n],[rj] it has 



been shown that the level statistics depends not only on the parameter r but also on the 
number of impurities Ni . This is clear in the superclean case. Then the density of states 
has narrow peaks centered around the unperturbed energy levels (|2.3|). The form of these 



peaks is Gaussian when the number of impurities is large ; however as shown in ||14| , |12f these 
peaks have a different form when there is only a single impurity. 

In the clean case u < 1/r < A, the level statistics also depends upon the number of 



impurity iVj. If iVj — > oo and if the random potential is a white noise, it was shown in |13 
that one can apply standard random matrix theory. However if the number of impurities 
TV, is smaller than a certain number N ic (N ic > (Ep/A) 1 ^ 2 ), Koulakov and Larkin [11] have 
found that the spectrum is of a different type. 

In the present article, we have defined the dirty and the clean limits by iV, > N ic , 
and Ni < N ic respectively. (These definitions differ slightly from the usual ones, quoted 
hereabove.) 

When there is only one impurity in the vortex core, with a radius £, the impurity is 
located at a position r^ ~ £. We are thus in a limit in which kpri » 1, since the correlation 
length is £ ~ h 2 kp/mA, and we have thus kp^ ~ k F /A ~ Ep/ A >> 1. Therefore, in the 
clean limit, this condition is always satisfied. In the dirty limit, the number of impurities 
increase and the typical value of the quantity kpti is reduced. 

In that limit the Bessel functions in ( |2.4| ) may be replaced by their asymptotic expansions 
and the matrix element A l nm , produced by the scattering site % between two states n and m, 
is given for kpti ^> 1 by 

A\ m ~ Ce i( - m ~ n ^ sm(2k F r t - ^±^ir) (2.5) 

where c = ^^(-M. 

The corresponding secular equation for the clean case was solved by Koulakov and Larkin 



TT| who have found that, after averaging over the location of the impurity, the density of 



states p(E) was given by 



p(E) = -sin 2 (— ) (2.6) 



(normalized to one after integration on the interval < E < uq). 

In the dirty limit, kpTi is not as large and it is necessary to take into account the sub- 
asymptotic behavior of the Bessel functions, and to sum over the impurities. Since the 
matrix elements A l nm have an oscillatory behavior, it is natural to regard them in the dirty 
limit as the elements of a random matrix. However it is necessary to keep some of the 
structure exhibited by the clean limit (|2.5|) in mind. 

In the dirty limit or in the moderate clean case, a random matrix theory has already 
been considered for the superconductor-normal interface, and for a quantum dot with a 
superconductor boundary. p!3|I5|. The excitation levels of the superconductor vortex have 
also been analysed in this way [[H|]. The random matrix ensemble, for such cases, has been 
suggested to be invariant by the symplectic group Sp(N). 

If the full Hamiltonian is treated as a random matrix in a symplectic ensemble Sp(n) , 
which is valid in the extreme dirty limit, the density of states has been conjectured |l5|JlT 
to be 



P{E) 



sin(^) 



"kE^ 



(2.7) 



This differs markedly from the clean result Q2.6|) . 

It is therefore challenging to find the crossover behavior between these two limits, and it 
might even be important for the transport problems associated with the excitations of quasi- 
particles inside the vortex core |14|. In order to achieve our goal, instead of treating the 



full Hamiltonian as random, we consider here only the scattering by impurities as a random 
matrix but not the full Ha miltonian. 

In order to motivate the use of the symplectic structure in this problem it is instructive to 
examine the form of the matrix A l nm . For a given impurity i, let A be the matrix whose 
elements are < n|74|m >= A l nm . The state index \m > varies over 2N levels. When N = 2, 
for example, the states \m > are | — 1 >, |0 >, |1 >, |2 >. Since E® = — uj n (n — i), we have 
E = -Ei = \uj q and E 2 = -E_ x = 
in the order |0 >, |2 >, |1 >, | - 1 >, 
2x2 blocks, we write 



-u {n- 2/ 
Labelling the lines and rows of the matrix A 
in order to split it into even-even, even-odd, odd-odd 



-fcuo. 



.4 



( <0|A|0> 

< 2|A|0> 

< 1|A|0> 
V< -1L4|0> 



<0|A|2> 

< 2\A\2> 

< l\A\2> 
< -1U|2> 



<0|A|1 > 

<2|A|1 > 

< 1\A\1 > 

< -1\A\1 > 



<0|A|-1 > ^ 

<2\A\ - 1 > 
< lU|-l > 
<-l|A|-l>/ 



(2.8) 



Since < 0L4|2 >= A^ 
have used the J- n (x) - 
N, the structure 



A 



2D 



< 2\A\0 >, and < 2\A\1 >=< 0\A\ - 1 >, etc. from 



we 



■l)™J n (x)), it is easy to see that the matrix A has, for arbitrary 



A 



a b 
b* -a T 



(2.9) 



where a is an N x N hermitian matrix a^ = a, and b is an N x N symmetric complex matrix 
b T = b {b T is the transpose of b). 



The number of degrees of freedom in the random matrix A is thus 2N 2 + N, ( N 2 for 
the hermitian matrix a and N 2 + N for the complex symmetric matrix b). The number of 
generators of the symplectic Sp(N) group is indeed also N(2N + 1). 

The symplectic structure of A is exhibited by the algebraic relation 



A T J +JA = 



where J is 



J 



1 
-1 



(2.10) 



(2.11) 



This structure implies that the eigenvalues are real and pairwise opposite, giving a chiral 
structure to the eigenvalues. The Lie algebra X is diagonalized by the symplectic group 
G G Sp(N). For example, for N = 2, we have 



/A, 



X = G ] 



A, 



■Ai 



V 



\ 



-a 2 ; 



G 



(2.12) 



If we mutiply the matrix A by % = ^/—T one recovers an element of the Lie algebra of 
the group Sp(N). If we consider now that the position (0j, r«) of the impurity is random, it 
is natural to take < n|A|m > in this Sp(N) invariant random ensemble. 

For instance when there is only one impurity (i = 1), the matrix A has a periodic 
structure, with alternating signs for each 2 by 2 block. For example, in the case of N = 4, 
we have for the full Hamiltonian H = E° + A, (E° is a diagonal matrix) 



.4 



/ s —s—c c \ 

-s s c —c 

-c c —s s 

V c — c s —s J 



(2.13) 



where s = C sin(2&;i7a) and c = C cosfikpa). ( In the single impurity case the phase can be 
set equal to zero by a choice of gauge condition ). This even-odd structure has been studied 
earlier in fll^llfl. The periodic structure becomes clearer if we write the matrix A = (A nm ), 
( f2.4| ) for N=2 with the approximation ( |2.5|) , 



.4 



( s c —s — c\ 

c —s —c s 

—s—cs c 

\—c s c —s J 



(2-14) 



in which rows and lines are in the order |2 >,|1 >,|0 >,| — 1 >. The determinant of this 
matrix is a kind of Toeplitz determinant (constant along parallel to the anti-diagonal). We 
have omitted the diagonal part E®. It is clear that the eigenvalues are periodic, depending on 
the odd-even parity of n. When we consider several impurities, a sum over impurities J2i A l nm 
should be taken. In the clean limit, we replace the quantities s and c by s = Yli C sm(2kFTi) 
and c = J2i C cos(2kpri). The random average over the positions of these impurities gives 
the density of state (g77|) |TTJ. 



III. ITZYKSON-ZUBER INTEGRAL AND ITS EXTENSION 

In the previous section, we have discussed how we are led to study the crossover from the 
clean case to the dirty case. We now divide the matrix SiJi^nm m ^° ^ wo parts. One part, 
denoted by H Q , corresponds to the clean case, namely the matrix has a periodic structure 
as in fl2.13|) , like a Toeplitz matrix, but with a sum over i. The other part is the remaining 
difference between the matrix Y,i=i A l nm and Hq, which is considered as a random matrix V. 
Thus we have 

H = H + V (3.1) 

The matrices H and V are both symplectic, represented by the symplectic group Lie algebra 
Sp(N), which has the form ( |2.9| ). 

In the analogous, but simpler, problem of a unitary invariant random perturbation, (H 
and V were then both complex Hermitian matrices), a technique to study the crossover 
from the non-stochastic Hq to a pure random matrix V has been developped earlier |2| -fl|. 
Here we have to consider the symplectic case. As in the unitary case, we have to integrate 
over the Lie group which diagonalizes the matrix H. In the unitary case, the first step, 
the integration over the unitary group, was introduced in the study of two coupled random 
matrices by Itzykson and Zuber. 0. We now summarize the formulae of integration over 
Lie Groups, which are useful for the subsequent discussion. We assume that the probability 
distribution of V is a Gaussian, and that Hq is a fixed non-random matrix. [It is possible to 
generalize it to some non-Gaussian distributions for V ; for instance in the appendix B, we 
have considered the case of a general hypergeometric function whose argument is a matrix.] 

The Itzykson-Zuber integral is an integral over the unitary group U(n), 



e 



Jr(uauH) du = C„ y e **') (3.2) 

lu(n) A(o)A(6) v ; 

where A (a) = I1j<;(oj — clj), and similarly A (b), are the Vandermonde determinants of the 
eigenvalues of the Hermitian matrices a and b. The constant is found to be Cn = IljLiO' - 1)- 
By expressing det(e ai6j ') as an alternating sum over the symmetric group S n , one may write 

e <^y a \ b>du = C N Ee ^\ 6 2T ( 3 - 3 ) 

u(n) A(o)A(6) v ; 

where Ad(u) ■ a = uau~ x is the adjoint action of u G U(n) on the matrix a. In this formula 
the matrix a is diagonal, a = diag(ai, ..., a n ), w is a permutation and e(w) is its signature. 
We have used the notation < a\b >= tr(ofo). 

This formula is a special case of a more general formula due to Harish-Chandra ||. Let 
G be a compact connected Lie group. Then Harish-Chandra's result reads ||6|, [T7| . |T8 



f p <Ad(gya\b> f]n _ T, w& w^w)e <w - ab> 
Jg 9 ~ A(o)A(6) [ ' 



where a, b, h belong to a Lie algebra h, and for any H G h, 
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A(H) = J] <*(H) (3.5) 

aeA + 

A + is the collection of the positive roots, w is the finite reflection group, called the Weyl (or 
Coxeter) group and h is the Cartan subalgebra. 

The Harish-Chandra formula has been interpreted more recently as an integration over 
an orbit, with a symplectic structure. This structure implies that the saddle point method 
is exact, provided one sums over all the critical points HJHj. This is also related to the 



localization theorems [20]. 



As an illustration of the semi-classical nature of this formula, we consider again the 
integral over the unitary goup U(N), 

1= f dge tr{a9b9t) (3.6) 

Representing an element g of U(N) as 

9 = 9oe lX (3.7) 

where X is an element of the Lie algebra, we expand the argument of the exponential, up 
to second order in X, and impose the stationarity condition, 

tv[ag iXbgl] + tr[ag b(-iX)gl] = (3.8) 

for any X, i.e. 

[b,glag ]=0 (3.9) 

If a and b are diagonal, it implies that go should be a unitary permutation matrix, g = p. 
Then, the integral / becomes 

I = ^2 e tr(ap ^ _1) I dX exp[ti(--a P X 2 bp- 1 - -apbX 2 p~ 1 + apXbXp)} (3.10) 

p J 2 2 

Performing the Gaussian integral over X one recovers the Itzykson-Zuber formula. (This 
is not a derivation of course, but a way of verifying that for this integral the one-loop 
approximation is exact provided one sums over all the saddle-points). 

A similar saddle point technique may be applied to the symplectic case. We take 
g G Sp(N) in fl3.6p and for a and b diagonal matrices with chiral eigenvalues, a = 
diag(ai, ..., a N , — oi, ..., — Ojv) and similarly for b. Then the same calculation yields 

I = c det[2sinh(2a I 5 i )] 

Hl<i<j<N( a i — a ])(H — b 2 ) l\i<k<N( a kbk) 

and this result will be repeatedly used below. 



IV. DENSITY OF STATE 



Using the method developped for the unitary case pj-[g,f|], we consider now the density 
of states with an external source matrix. Indeed for a random Gaussian V, the resulting 
probability for H is a Gaussian with an external matrix source A linearly coupled to V. We 
assume that the matrix A belongs to the Lie algebra of the symplectic group and, without 
a loss of generality, we may take it as a diagonal matrix, A = diag(ai, ..., ojv, —a-i, ■■■, —o>n)- 
We take this A as the unperturbed H Q . The probability distribution of the random matrix 
V is assumed to be Gaussian. 

The density of states p(X) is given by 

p(A) = -^<f>(A-A a )> (4.1) 

The probability distribution P(M) of the random matrix M is 

P(M) = exp[-trM 2 - tiMA] 

= exp[-£A 2 -tr(< ? A< ? U)] (4.2) 

where A = diag(Ai, ..., Ajv, — Ai, ..., — Xn) and g is an element of the symplectic group. (Note 
that we use here a normalization which differs from our previous treatment [0,|]|| of the 
unitary case). With the present normalization the edge of the density of states becomes of 
order yN ; in the large N limit, the support of the density of states lies in an interval of 
order [-y/2N,y/W\. 

The density of states is obtained as the Fourier transform of the evolution operator 

/OO ri+ 
.. H e ^(*) (4-3) 



and using the symplectic measure A 2 (A) with A(A) = Il(^ 2 — A 2 ) II Xk |2l| , and the previous 
Harish- Chandra integral formula ( |3.11| ), we obtain 

iV a=l J0 j=l i.ll<i<j<N\ a i ~ a j) li.l<k<N a k 

x det[sh(2A i a j )] (4.4) 

(up to a factor 2 in (|3.11|) absorbed in the coefficient C). 



By the reflexion symmetry A» — > — A«, and by the exchange symmetry between Aj and 
Aj, we can extend the integrations over the Aj's from — oo to oo . The change of sign for 
itXi can be absorbed since we can change t — > — t by parity. Then the expression for UA{t) 
simplifies to 

— V f°° TTHA n^!<^^j ~ A i)ni<fc<7V Afc y <x?+it\ a +2Y.*' a,X,. 



«*W = ^E/ n^ n ,Z I "'^"V ^.+''^<^ (4.5) 

JV a =l J -° i=l i.ll<i<j<N\ a i ~ a j) lll<k<N a k 

Before proceeding to the explicit integrations over the Aj's in U^t), let us investigate 
the expression for the density of states p(X) , in the absence of any external source, which 
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one may then calculate by the usual method with orthogonal polynomials . The N-point 
level distribution is given 



p N (X 1 ,...,X N ) = C J] (^-A?) 2 II ^exph^A 2 

l<i<j<N l<k<N 



(4.6) 



Using orthogonal polynomials method ||10|1 , this distribution may be written as 

pjv(Ai, ..., X N ) = det[K N (X h Xj)] (4.7) 

with the kernel Kn(\{, Xj) given by a sum of Laguerre polynomials, and in particular the 
density of states p(X) = Kn(X, X) is equal to 



1 ^ n! 



™ = vXrifaj#\*)#\w-»- 



(4.8) 



n=0 *- \2 



Ln ( x ) are associated (generalized) Laguerre polynomials, orthogonal with the normaliza- 
tion 



6 X 2 Lin (X)L/rn (XjCLX — : ^nm\ 

n\ 



(4.9) 



leading to L { o /2) (x) = l,L^ /2 \x) = 3/2 - x, L^ 2 \x) = E?= (-l) r [r(n + 3/2)/T(n -r + 
l)r(r + 3/2)}x r /r\. 



The asmptotic behavior of these Laguerre polynomials in the limit N —> oo is [p3|1 



e3* sin(2 VWx) + 0' 



4V) 



With the Christoffel-Darboux identity, we have 

N-\ 



N 



E 



„=o r (f + ™ 



^—^I^\x)I^\y) 



r(| + iv) x-i/ 

and using the large N aymptotic behavior ( |4.1U| ), we obtain 



re 



N-l 



n= ni + n 



f-—^L ( n 2 \r)L { n 2 \r)} 



xe-^LfW^Ur) - Lf^x^Lf {x)\ 



l r VN r 

sin 



7T \/X 



N' A-s/Nx 

~ -[1 == sin(4v / iVx)l 

tt l 4\/JVa; 



sin(4\/iVx)] 



(4.10) 



(4.11) 



(4.12) 



Thus, putting x = A 2 , we get 

vr Ay/NX 

This result is identical to (|2.7|) , when we substitute o; = ^ and £? = X/y/N. 

We now return to the integral ( }4.5| ) over the Aj's in the presence of the external source. 
We first replace t — ► It and following @|§, for fixed a, we denote by 6; the sum hi = ai + itSi t0l 
. One then uses the integral 

/ dA.niA^-A^nv-^^^n^-^n^ (4.14) 

(after the translation Aj — > Aj + fej the result follows easily from antisymmetry under permu- 
tation of the fej's, parity in bi, and counting the degree of the resulting polynomial in those 
6's). 

Using the normalization Ua{0) = 1, and writing the sum over a as a contour integral 
around the a 2 's in the complex w-plane, we obtain 

1 / du £ ((Vi + if) 2 -ag 1 it.^ +2a ^ 



U *® = jf> 2S II ( U -aj) (^ + ^-^ + 7i )e " ' (415) 

Fourier transforming U A (t) we obtain the density of states p(A) 

p(A) = £u A {t)e-*»^ (4.16) 



Let us verify the consistency with the orthogonal polynomial result (|4.8| ) when the external 
source vanishes. Then 

Uo{t ) = U ^ + ^ 2N \ (1 + « )«-*»** (4.17) 



and 



dt -- 2it\ 



/oo w/- 



i(-l) N r°° dx r du , x , »r,i 1 



27V+1 x -a 2 -2ia;A-M+2AVM /^ io\ 



2 



./V J-oo 27T J 27U y^ W + X 

(we have shifted ttot = x + i^/u). The integration contour is a small loop around the origin 
u — 0. Expanding l/(w + x 2 ) = x~ 2 (l — -^ + • • •), and noting that 

— —L— e - u+2X ^ = ^\ (-l) n \L$\\ 2 ) (4.19) 

27iiu n + 1 ^ T(| + n) v ; v ; v ; 

and 
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~° dx^ 2n+1 2_ 2ixX _ in\X (|) 2 X 2 



OC 



x j n +L e -x"-2 lX A = lLl£>(\*) e -*' (4.20) 

2ir 2y/n 



we v erify the consistency of the integral representation (|4.18|) with the orthogonal polyno- 



mial result ( |4.8| ). The integral representation ( 4.20J ), gives an easy way to recover the large 
n asymptotic behavior ( (4.11 ) of the Laguerre polynomials through the steepest descent 



method. Note that the chiral invariance leads as expected to an even density of states. 
In the presence of the external source, the density of states becomes 



N J-oo 2n J 2m y/u ±r x u — a^ u + x z 

(We shall analyze this density of states as the limit of A — > \x of the kernel Kn(\, fj) in the 
next section). We note here that if we make the change of variable, u = v 2 , we have 

p(x) - t{ - l)N ' 1 r - j dv f\( x2+a * ) x c -*-**x-*+2xv f422) 

P{X) ~ N y-oo27rJ2vr S 1 = V^_ G 2^ 2 + x2 e ^ ZZ > 

where the contour encloses all the ±o 7 . For example when N = 1, we obtain 



... A sinh(2Aai) n i X 2 

pW = irr — e ' 



A e-( A -^) 2 _ e -(A+a 1 ) 2 ] (4>23) 



4^01 



Near the origin A = 0, the density of states p(A) becomes p(A) ~ — exp(— a 2 — A 2 ). Thus 
the density of states is an even function and it vanishes near the origin as A 2 . 



V. CORRELATON FUNCTION 



The n-level correlation functions are expressed as the determinant of the kernel K^(X, (J>), 
exactly as for Hermitian random matrices in a source. Indeed the formula of ( |4.14[ ) is quite 
similar to (3.8) of ||. Therefore, all the derivations of the Hermitian matrix model in |§f|||, 
may be repeated almost identically for the Sp(N) case. 

For instance, the two-level correlation function p^ 2 '(X, p) is obtained as the double Fourier 
transform of 6^4(^1,^2), 

P {2 \Kp)=JJ^e-^-^U A (tut 2 ) (5.1) 



in which C/a^i,^) is defined as 



I N I 



U A {hM) =< T7 E ^T7 E ****** > (5-2) 

iV a=l iV (3=1 



11 



Using the integral formulae ( |3.10| ) and (|4.14 ), we find 



TT u +, 1 ^ {a a + ih) 2 - (a + it 2 ) 2 (a a + itif - a* 
UA[t u t 2 ) - -= ^ , 2 _ 2) 11 -j— 2 

x -j-j ( a /3+^i) 2 ~ a 7 (1 , it^ | ^ 2)c „ t 2„ f 2 +2nMi+2 ^ f2 



(5.3) 



(wa have subtracted the term a = (3 in £/a(£i, £2))- The summation over a and /3 is expressed 
by contour integration over two complex variables u and v, 

u A (t 1} t 2 ) = — f ^-^ n — 7^zr^\ — - n 

x 



N 2 I (27Ti) 2 ^ (« - a 2 ) 11 ,, - a 2 

1 1 iti s ,_ ito 



a + —)(i + —) 

{u - v) [{y/u + ihf - (y^J + it 2 ) 2 ] 

[{^+it l y-u]\{^ + it 2 y-v} 



x (W - V) [Wu + iti) 2 - (y/v + Jt 2 ) 2 ] r - t 2- t 2 +2l ^lt 1+2 i^t 2 



We now shift ti — ► ti + «a/u and t 2 — > t 2 + «>/^ n °t e that 

(«-t;)pi) 2 -(^) 2 ] [(i<i) 2 -u][(»fa) 2 -t;] 

[(ztO 2 - * 2 ) 2 - «] " [(^i) 2 - v}\{it 2 y -u}~ 



(5.4) 



(5.5) 



and divide UA(ti,t 2 ) into two terms. The first one, which comes from (-1) gives the product 
of the two density of states p(A)p(/i), and the second one is a product of two integrals over 
ti,u and t 2 ,v respectively. Therefore, we end up with 

p (2) (A, n) = K N {\, \)K N {n, p) - K N {\, n)K N (ji, A) (5.6) 

with 

x (^ ); (5 - 7) 



(the contour in the w-plane encircles all the a 2 ) . From this expression, we recover the 
previous result ( f4.15| , [4.16l) for the density of states as p(A) = K N (\, A). 



The n-point correlation can be analyzed as in |J and leads here also to a determinant 
form, 

p(\ 1 ,---,\ n ) = det[K N (\ l ,\ j )} (5.8) 

for an arbitrary external source. 
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We now discuss Dyson's short- distance universality, within our model, i.e. when we vary 
the external source A, through the explicit expression of the kernel K N (\,fi). 

We return first to the sourceless case,a 7 = 0, and consider the large N limit of K N (\,ft) 
for A, /i are order of one. We scale t — > yNt, and u — > Nu. The kernel is then 



K N (\,p) = --L fp: /^. e -^C/(t)-H»W)( « ) 1 (5.9) 

v ' p; y/NJ2mJ 2vr K x/u J t 2 + u v ; 



where f(t) = t 2 + 2ifj,/yN — 2 In t and g(u) — u — 2\\fuj yN + In w. In the large N-limit we 
have the saddle points, t { p = 1 - ifi/2\/N, tf> = -1 - iy,/2y/N,v$ = -1 + i\/2y/N and 
u^ 2 -* = — 1 — ifj,/2y/N. Adding these saddle-points contributions, we obtain 



V ' w tt [ 2y/N(X-fi) 2v/iV(A + ^) J V ; 

When A — > /i, we recover the expression (^4.13|) of the density of states. 

In the presence of the external source, we find a similar expression to Eq.( |5.9| ), with 



/(*) = t 



2itfi 1 a 2 






, , 2\ x [u 1 ^-^ . a 2 . 

v/JV N ^ N 



Using the definition , f'(t c ) = 0,g'(u c ) = 0, of the saddle-points we find 



= IMW^A-a*)) _ sin(W*(A + „)) 
V ^ tt [ 2VN(\-fj,) 2%/iV(A + /x) J V ; 



where t c o is a solution of 



1 



-V - = 1 (5.14) 

7=1 l c0 + AT 

Here we have assumed that the order of a 7 is a 2 < O(N). The derivation of ( |5.13| ) is the 
same as ( |5. 10| ) . The difference is just a normalization, due to the change of the saddle-point 
t C Q. Note that the support of the density of states is inside the interval of —\JN and \JN . 
Therefore, universality holds provided that the eigenvalues of the external source matrix are 
located in an interval of the same order of magnitude ; ( it has to be of same order as the 
support of the density of states for the zero external source case). In this universal regime for 
the correlation function, the energies A and /x are assumed of order one, namely X/\/N — ► 
in the large N limit. Since A can be order of y/N, the universal behavior of ( |5.13| ) appears 
only near the origin, as exepcted, and is differnt from the non-universal bulk behaviour. 
The universality found here, is similar to that found in other chiral random matrix models 
studied in ||; there it was the density of states which took a universal form near the origin, 
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and became independent of the external source or of the non-Gaussian distributions (as seen 
by a rescaling of the energy) . 

When the a 2 spread over an interval larger than N, the universal form of ( |5.13| ) does not 
hold any more, since the saddle-point method of (|5.14j) breaks down. In our formulation, 
the amplitude of the random matrix is fixed. We now change the strength of the external 
source. For this purpose, we introduce a parameter C, which determines the strength of the 
external source Hq, respective to the random potential V. Let us consider the example of 
the external source, 



Ctan 



(27-1)71-, 



L 2(2iV+l) J 
where 7 = 1, 2, ..., N, and C is a parameter. When r y/{2N + 1) << 1, we have 

a = ° n (7--) 
7 (2iV + l) 17 2 j 



(5.15) 



(5.16) 



and the eigenvalues of the external source are equally distributed, and thus reproduce the 
excitation spectrum inside the clean superconductor vortex, if we identify Tiujq = Cn/(2N+1) 
in (|2.3|) . Since we are interested in the behavior near the origin for the energies A and \x in 
fGv(A,/i), this linear approximation of the tan(x) ~ x is valid. 
More generally, using the formula, 



SL . n 2r (2m — 1)7T 1N 

n(x 2 + tan 2 [ ^ oAr ; : ' ]) 



1 



m=l 



-2(2iV + 



2(2iV+ 1) 



[(1 



x) 2N+l + {l _ X2AT+11 



(5.17) 



we may replace the product involving the a 7 in 



and obtain 



™">~*/s/£ 



;i + {,f+ l 



±\2N+1 
C> 



W2i\2N+l 
C I 



'I _ V^2JV+1 



1 



it 



t 2 + U y/u 



)e 



-t 2 -2it/i-u+2 y / uX 



(5.18) 



This expression is exact for finite N when the a 7 are given by ( p,15| ). It is easily seen 
that when C — > 0, we recover the previous external source free result ( |5.9| ). 

For non-zero C we have two different cases according to the size of C, 
i) if C is of order \/N, we recover the universal behavior of ( |5.13| ). (We simply scale t — > \/Nt, 
and u — > Nu, and use the saddle point method). We recover then the result ( |5.13| ), with 



t 



cf) 




c 



N 



(5.19) 



ii) in the second case C ~ O(N), we change u = v 2 in ( |5.18| ) : 



*«*■") -4/s/5 



;i + i) 2N+1 + (1 - i) 2N+1 



'l _|_ vv\2N+l 



'I — vv\2N+l 



X 



it 

t 2 + v 



-t 2 -2itii-v 2 +2v\ 



(5.20) 
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In the large N-limit (|5.20|) reduces to 



^"> = -*/£/ 



1 r dv r dt cosh(^) it 



2tt cos(^) £ 2 + v< 



If we rescale t — » -y/iVt and t> — >• v^Vv poles appears in the u -plane at v — C(m — |)7r/2iV 3//2 
for which cos(2N 3 ' 2 v/C) = ; m is an integer. 

Within the residues of these poles, we have the Gaussian factors 

JVC, A 12 ,C(m— i)ir , 



There are thus two different regions for A. We first consider A ~ 0(y/N), i.e. the bulk 
case. Then X/y/N » y/N/C, and the saddle point of t is t = —iX/y/N. We set A = /i, 
and for this this saddle point, we have the density of states p(X) = K N (X, A) 

Acos(- 27VM 



2y/rrN m _\2 + ( c ( m ^ 



2N 

x e -(^iv^- A ) 2 (5.23) 

Note that the zeros of the denominator are cancelled by zeros of the numerator 
cos(2NX/C). When C/N » 1, the exponential factor of (|5.23|) damps the result and 
we may set A = C(m - l/2)n/(2N), then 

A " S( ^> ->-!)-£ (5.24, 



-A 2 + (^#^) 2 ~ C 

Thus we find 

P(A) = 7^7=E^ ( ^^ A) (5.25) 

which is a sum of Gaussians : the density of states is just sum of Gaussian peaks around the 
levels A = C(m — 1/2)tc/2N. This result is valid provided A is not too close to the origin. 

When A is near the origin, namely A ~ 0(1), we have X/y/N ~ y/N/C since C is order 
N. In this case, we have to take into account the term cosh(2iV£/C) in the saddle-point 
equation. The saddle-point becomes to = —iX/y/N ± y/N/C. Using this value in ( |5.21| ), we 
find 

A -4- i^- 

"W = ^B-in _ (A + i g )2+ + ( | (n _ |K) > -'""-^' (5-26, 

Near A = 0, the density of states behaves as p(A) ~ A 2 . For example, when C = N, we have 
p(X) ~ 7rA 2 exp(-vr 2 /16)/(l + vr 2 /16), and p(A) has peaks at A = ±vr/4, ±3vr/4, .... Thus, 
when C ~ N, which means that the random potential is weak compared to the external 
source, the density of states takes a form which is different from the universal form ( |5.13j ). 
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VI. CROSSOVER FROM THE CLEAN LIMIT TO THE DIRTY LIMIT 

In the clean case, when k F ri ~ k F t; » 1, we approximate A % nm by (|2.12|). Then, the 



periodic spectrum is obtained 

E n = -hcu (n-^ + 2(-l) n ) (6.1) 

When the number of impurities inside the vortex N t ~ 1, we are in the superclean case, and 
z in ( |6.1|) is a function of the position of this impurity. If 1 < Ni < iVj c , we are in a clean 
case, and we need to average over the different values of z. 

In the previous section, we have discussed the case where the deterministic term E n = 
—Tiuj{n — 2) is coupled to the random matrix A. Here we consider for the external source 
the matrix whose eigenvalues are ( |6.1| ). 

We may then use the formula 

W + tan2( w^T)>] = i(iW l(1 + x)2lN+1) - (1 - 1)2, " +1,| • (6 ' 2) 

which is similar to ( |5.17|) . ( We assume that N is odd here). We divide this product into 
two parts, r-odd and r-even . The r-even part is obtained immediately from ( |6.2| ) as 

Ar - 1 rvr „ 1 



TT [ x 2 + tan 2 ( ' - )] = - N [(1 + x) N+1 - (1 - x) N+1 ] (6.3) 

Dividing (|6.2| ) by (|6.3| ), we obtain the expression for the r-odd part, 

fiy + tan2 (^ry )] = \iv + *) N+1 + a - ^ +1 ] ( 6 - 4 ) 



r=odd 



We now consider the eigenvalues ( |6.1| ) as an external source 



a 7 = [Cn/2(N + 1)]( 7 - I + z(-l)?), (6.5) 

and introduce z = Ctt/A(N + 1), z = —Ciiz/2(N + 1). Using the expressions ( |6.3| ) and 
TM, we write 



,X — IZq — iz %-} X — IZq — IZ 2 2/ r7T 



^_ )2 + tan2 , 



C r=even & 2(iV + 



x TT U X ^ + ^ ) 2 + tan 2 ( ,""* - )] 



" 4(iV+l) K C ' ~ [ C ' J 

x [( i + x-izo + iz )N+1 + (1 + x-izo + iz )N+1] (6 6) 
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In the large N limit, when m/(2(N + 1)) << 1, the left hand side of ( |6.5|) vanishes at 
x = -i(±r - 1/2 + z(-l) r )Cn/[2(N + 1)] . Noting that 

Y[(t 2 + a*) = IJ[* - «(*b + z)+ i „ fA j^^ ][t ~ *(zo + z)-i- t ' ' : " ' 



2(iV + l) JL v u ' 2(iV + l) J 
II[(t - i(z + z)f + C 2 tan 2 (^^)] (6.7) 



we have from (|5V7l ) (and the change of variable y/u = v), 



x 



X 



N J 2iriJ 2n t 2 + v 2 

iHt- Zo -z)r-(i-i( 



[(1 + U-it -z Q - z)) N+1 - (1 - U-it -z - z)) N+1 ] 



[(1 + £(« -z - z)) N ^ - (1 - i(v - z - ^)) 7V+1 ] 

[(1 + i(-g - * + z))^ 1 + (1 - i(-g - z + z)) N+1 ] 

[(1 + ±(v - z + *))" +1 + (1 - £(v - ^o + ^)) 7V+1 ] 



(6.8) 



As in the previous section, we have to distinguish two different cases, i) C ~ O(N) and ii) 
C ~ C^-y/iV). There are also two regions A ~ 0{yN) and A ~ 0(1). 
When C ~ O(N), we exponentiate the factor in the bracket of 



sin[^(— it - z - z)]cos[^(—it - z + z)\ 

sin[^(w — z — z)]cos[^(v — z + z)] 
sm[2g-(-it - z )] - sin(^z) 
sm[^(v-z )]-sm{^z) 

Since z = Ctt/[4(A^ + 1)], we have 

1 r dv f dt cosh(^) + sin(^z) it 



(6.9) 



1 r dv r at cosh{^) + sm{^hz) it 



When z = 0, we obtain the same expression as in the previous section ( |5.21| ). 



Poles in the integral over v are present at v = z$ + z + Cnn/N and v = Zq — z + (2n + 
l)(Cn/2N) from (|6.8|). The residues i? for the first poles are 



p _ /■ 1 r?Y Sill [f (~^ ~ ^0 ~ gjWgH^ ~ *0 + Z)\ rfi11 s 

1 j [t 2 + (z + ^+^) 2 ]cos[f(^ + ^)] Ib - Uj 



When A = /i ~ 0(VN) } the saddle point becomes to = —iX/vN and we set this value in i?. 
The second pole gives a similar expression. Furthermore, if C/N » 1, we can approximate 
A ~ zq + z + Cnn/N. Then, we find that the residue R in ( |6.10| ) becomes R = N/2C. 
Therefore, we have the density of states for C/N » 1, and A ~ 0(y/N), 
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P(A) 



N 



AC^ 



El 



-( Z0+Z+ 2™L- 



-xf 



+ e 



-(zo- 



C(2n + 1)7T 

2iV 



-A) 2 



(6.12) 



which is a sum of Gaussian distributions. 

In the clean case we have several impurities and one may average this density of states 
over z . As discussed in [11], the measure for this ^-average is 



po(z) = Acos 2 [ 



2Nz 
~C~* 



(6.13) 



where A = AN/Cir. The matrix elements are represented by a quaternion, and the measure 
is isomorphic to the uniform measure over the four-dimensional sphere s 3 , Sp(l) ~ S 3 . In 
spherical coordinate this leads to ( |6.13| ). When C/N » 1, the integration of ( |6.12|) over z 
with the measure ( |6.14| ) yields 



< p{z) > 



r e-^+'+^-V Acos\ 2 4^)dz 

J — oo 



c 



jix iN 2 AN 

4^[l-e-^cos(— A)] 



(6.14) 



in 



which we have used N/C « 1). When C/N — > oo, the exponential factor 
exp(— AN 2 /C 2 ) becomes one and the density of states becomes sin 2 (2iVA/C). This is 



[Till . When C/N is finite, at the energy 



indeed the result for the clean case found in 
A = mrC/2N = nhuo, the density of state becomes non- vanishing, except at the origin. At 
the origin A = 0, we have a zero due to the factor it = A in (|6.11| ). It may be interesting to 
note this deviation from zero has been observed in the numerical work of [ |ll|] . 

When C becomes order of y/N, we may apply the saddle point method in (|6.8|). We take 



the saddle point equation (|5.14j ) for a 7 = ^{n — | + (— 1) 7 5). However, the second and 
third term (—1/2 + (— 1) 7 5) can be neglected since they remain of order one, compared to 
n which can be order of N. Thus we may approximate a — 7 by a 7 = 7r 2 C 2 n 2 /4iV 3 , which 
is independent of z. Therefore the average over z does not change the result for the density 
of states, and we have a universal result for the kernel and for the density of states ( |5.13| ). 
It is possible to take the average (|6.10|) by the formula, 



I ^ a + smz ^^ = 1 _ _ ^3T)i 

■5 6 + sinz L 2 v /v n 



(6.15) 



Then, changing the contour integration in the f-plane to the integration near the saddle- 
points for v in the large N limit, we have 



K N (\,/j) 



CAtc r dt r dv 



t 



x sin 



AN 2 J 2tt J 2ttH 2 + v 2 

'2NV . - t 2-2itn-v 2 +2vX 



2Nt 2Nv 

)[(cosh(— ^)-cos(— — )) 



C 



C 



C 



(6.16) 



The term sin(2Nv/C) comes from the sum of integration paths in opposite directions. We 



assume here that A and /i are of order one, and 0(A) ~ 0{-^). After the shift t — > \/Nt, the 



saddle-point is at to = —i^/VN± vN/C and v = X/yN ± iyN jC for the terms in ( |6. 16|) 
Inserting the values of these saddle-points, we have 

CX 2 , _***' ANX S 



^[1 + e ^cos(^^)] 



4iV(A 2 + -^) L y C 

+ i^T|l) e ^ sm( — > (6 - 17) 

The above expression reduces to that of the dirty case when C/N — > 0, since the first 
term is then dominant. Thus ( |6.17D gives the correction to the dirty limit, and it is valid for 



C ~ N. 

In the large C, C » N, it is easy to recover the clean case from the expression ( |6 . 1 6| ) . 



The saddle-point is given by t = —ifi and v = X. Putting these values in ( |6.16[ ), and taking 
fi — > A, we obtain immediately 

p(X) = — sin 2 (— ) (6.18) 



Thus, ( |6.16| ) gives the result both for the dirty case and for the clean case when one 



varies the parameter C. 



VII. DISCUSSION 



We have discussed a random matrix theory for the energy levels inside a superconductor 
vortex and investigated the crossover from the clean case to the dirty case. 

The technique involves a transposition of earlier results for the Hermitian case to the 
symplectic group Sp(N). An extension of the Itzykson-Zuber integral is presented (in 
an appendix). In the symplectic case we have to consider an external source matrix, 
which describes the spectrum of energies in the absence of impurities, whose eigenval- 
ues are of the form a = (a\, ..., ajv, — Oi, ..., — ajv)- In the symplectic case the measure 
A(ci) = n(a 2 — a 2 )]\ai replaces the Vandermonde determinant of the unitary case. In- 
deed, in the superconductor vortex, the eigenvalues appear in opposite pairs, i.e. with an 
exact particle-hole symmetry, as seen in the Andreev reflexion. We have then discussed gen- 
eralized integral formulae related to the Sp(N) group. In an appendix we have investigated 
hypergeometric functions of matrix argument both for the unitary and symplectic groups. 
This result may be useful for non-Gaussian random distributions. 

This random matrix theory is phenomenological. One phenomenological parameter, 
which can be viewed as the intensity of the disorder, is sufficient to describe the distribution 
functions of the random matrices. We have fixed this disorder parameter to one, which is 
why it does not appear explicitly in ( |4.2|) . Instead for the external source matrix A, we have 
introduced a parameter C in ( |5.16| ). We find the a universal formula for the density of states 
and the correlation functions for the different regimes of this parameter: a clean case C > N, 
a dirty case C < N, and a crossover region C ~ iV '. Nevertherless the identification of this 
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parameter C with microscopic parameters, such as the strength of the impurity potential Vi 
and the number of impurities iVj in (|2.4|) , requires a microscopic study which is outside the 
scope of this work. 
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Appendix A: Harish-Chandra formula for the unitary, orthogonal and sym- 
plectic group 



Let G be a compact connected Lie group. The Harish-Chandra integral formula is 

0JH 



a]b> , = E w ew(detw)e< w - a \ b > 
y ~ A(o)A(6) 



,<Ad(g)-a\l 



where a, b G h; h is a Lie algebra, and for any H G h, 

A(H) = 1] <*(H) 



(A.l) 



(A.2) 



A + is the collection of positive roots, W is the finite reflection group, called the Weyl (or 
Coxeter) group and h is the Cartan subalgebra of the Lie algebra of the group. In the general 
classification theory, the irreducible finite refelection groups are categorized as belonging to 
various types, A n , B n , C n , D n , ...., which are associated with certain compact Lie groups. In 
( p.4|) , det(u>) is simply ±1, since each w G W is an orthogonal transformation. 

Let us illustrate this result in the simplest case of the unitary group. The Lie algebra 
of £77 (n) , u = su(n), consists n x n complex skew-Hermitian traceless matrices. The 
complexification of u is g = sl(n,C), the Lie algebra of all n x n complex matrices with 
zero trace. The Cartan subalgebra h of g consists of diagonal n x n complex matrices 
H = diagi^hx, ..., h n ) such that hi + .... + h n = 0. We define the linear functional Cj on h 
by Cj(H) = hj, and the n x n matrix Ej^ which consists of 1 in the (j,k)th position and 
elsewhere. The linear functional a — e j — e^, j ^ k is a root of g with respect to h, i.e. 
A = [ej — efe : 1 < j ^ k < n]. One then verifies that ( |A.1| ) reduces the Itzykson-Zuber 
formula ( p.4|) in the unitary case. 

The orthogonal group O(N) consists of real N x N matrices u such that uu l = 1. The 
Harish-Chandra formula applies to compact connected groups G and we thus restrict our- 
selves to the special orthogonal subgroup SO(N), of orthogonal matrices with determinant 
one. The Lie algebra g = so(N) of G = SO(N) consists of all iV x iV real skew-symmetric 
traceless matrices. We need to consider separately the even case N = 2n and the odd one 
N = 2n + l. 

For SO(2n), the Cartan subalgebra h is the complex Lie algebra of all 2n x 2n complex 



block-diagonal matrices of the form p2 



H 



( ° 


h 






-h 
















h 2 






-h 2 






V 



o 



(A.3) 



J 



This matrix H is written as the direct sum of v's defined as 
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-1 J) 



H = hxv © h 2 v © • • • © h n v (A.5) 

Let [ei, ...,e n ] be the standard basis for i? n . A root system for SO(2n) is A = [±ej ± e^ : 
1 < j < & < n], which is a root system of type D n . For a = Cj ± e& G A + and the Cartan 
subalgebra H, we have a(H) = hj ± h^ and 

*W = II (h'-hl). (A.6) 

A fundamental Weyl chamber is defined as S = [(hi, ..., h n ) E R n : hi > ■ ■ ■ > h n -i > \h n \}. 
The group G(n) of permutations w of the set [— n, ..., — 1, 1, ...,n] restricted to w(—j) = - 
w(j), acts on the set of [h- n , ..., h-i, hi, ..., h n ] as 

w ■ (hi, ..., h n ) = (h w (i), ..., ^ W (n)) ( A -7) 

where we denote h-j = —hj,j = 1, ..., n. The Weyl group W consists of W — [permutations 
and even number of sign changes of [e\, ..., e n ]]; thus \W\ = 2 n ~ l n\. For H, we have 

w ■ H = h w( i)V © h w (2)V © ■ ■ ■ © h w (n)V (A. 8) 



Therefore the Harish- Chandra formula ( |A.1[ ) gives for a = a x v © o^f © ■ • ■ © a n v and 6 = 

b t v © M © • ■ ■ © &n« 

.t^r 1 ^ _ c *n, . E ^ G W (detw) exp(2 E ? =1 ^l^j f A q) 



SO(2n) '' , ' ! "" ni<j<fc<nK 2_a fc)( 6 J ~ 6 



with Csgh = (2n - 1)! I^ l (2j - 1)!. 
For SO(2n + 1), the matrix if is 

H = h 1 v®h 2 v@---@h n v@0 (A.10) 

and the root system for SO{2n + 1) is 

A = [±ej ± e k : 1 < j < k < n] U [±e i : 1 < j < n] (A.ll) 

Thus for a 6 A + and H, we have 

a(if ) = hj ± /ifc i/ a = tj ± e^ 

ft,j z'/ a = ej (A. 12) 

and 

n 

v(h)= n (^ 2 - ^D n ^ (A.i3) 

l<j'<fc<n j=l 

22 



The weyl group of SO(2n + 1) is W =[ permutations and sign changes of [ei, ..., e n ]]. \W\ 
2 n n\. 

The action of G(n) on the Lie algebra h is 



W ■ H = fl w (i)V © /l w (2)U © • • • © /i«,(n)^ © 

The Harish-Chandra formula (|3.4|) becomes for a = a\v@- ■ -@a n v ©0, 6 = b\V ( 



(A. 14) 

■©6„w©0, 



/■ e^- 1 ^ _ <?„ , Eu, eG( n)(det W )exp(2E7 = iw(a J -)& 3 -) fA , 

Ao^n ^ G(n) n,™^? - aim - m m_, a ,-&,- l J 



'50(271+1) 



ni<j<fc<n( a j ~~ a k)\Pj ~~ "k) ilj=i Oj'Oj 



where C G(n) = n?=i(2j - l)UJ^jl 

For the case of the symplectic group Sp(N) , we have 

fhx 



H 



h, 



-h 



(A.16) 



V -K J 

and ej(H) — hj. A root system for Sp(N) is 

A = [±ej ± e fc : 1 < j < k < n] U [±2e,- : 1 < j < n] 



(A.17) 



The Weyl group for the Sp(N) algebra is W = [permutations and sign changes of [ei, ..., e n ]]. 
\W\ = 2 n n\. 



Appendix B: Generalization of the Itzykson-Zuber formula 

The Itzykson-Zuber formula for Hermitian matrices a and b is again 

det(e aib J) 



e .r(gag 1 b) rf 



,tr 



U(n) 



y 



A(o)A(6) 



(B.l) 



where Oj and foj are eigenvalues of a and b, respectively. We wish now to consider general- 
izations of this type of integrals for which 



u(n) v ; A(o)A(6) 



(B.2) 



where ip n is a function of a matrix argument, and / is a real function. 

The Itzykson-Zuber formula of QB.1|) corresponds to /(a;y) = exp[xy] and ^ n = 
exp[tr(#a# -1 &)]. 

We may take for f(xy) an hypergeometric function p J r q (ai, ..., a p ; fa, ...,(3 q ; xy) . Note 
that f(xy) = exp[xy] = J r o(xy). Then , the corresponding function ip n {t) is also an hyperge- 
ometric function of the matrix argument p F q (ai+n—l, ..., a p +n— 1; /3j+n— 1, ...,/3 g +n— 1; £). 
For example, in the p = 1, q — case, we have 



23 



/ det(l - agbg- l )- a - n+l dg = Ma; a, b) 

JU(n) 



U{n) 

Hptf, TJrva,hX\ 

B.3) 



det(iJF (a; dibj)) 



A(o)A(6) 

where iJ-"o(a; dibj) = (1 — dibj)~ a . This is easily checked directly for the n = 2, g = U(2) 
case. We represent g as 

/ cos 4>e %ei sin (fie 102 \ m a\ 

9 ~ {- sin <j)e i9 s cos fa-Wi-**-**) J ^-4) 

with the measure J = cos sin 0<i0 Flf =1 d6i. 

More generally, the formula of the integration over the U(N) group may be written as 

p F q (a 1 , ..., a p ; j3 u ..., j3 q \ a, b) 

= / p F q (a 1 ,...,a p ;f3 1 ,...,flq;dgbg~ 1 )dg 

JU(n) 

_ det{pJ r q (ai - n + 1, ...,a p - n + l;/3i - n + 1, ...,(3 q — n+ I; dibj)) 
~ A(o)A(6) 

(B.5) 

This derivation of ( B.5|) proceeds by induction [p5fl . First, we introduce the zonal poly- 
nomials Z m (A) |24| which are homogeneous polynomials of degree m, which are symmetric 
functions of the n eigenvalues of the matrix A. From their definition they have the simple 
property that 

Z m ( 9 Ag- 1 B)dg = Z ^f Z ^ (B.6) 

in which the integral runs over the elements of a compact Lie group G with a Haar measure 
normalized to one. The coefficients of these polynomials are group- dependent, but they 
can be constructed from this property inductively. The polynomials are thus expressed as 
decompositions of products of Tr(A™)Tr(A™) ■ ■ -Tr(A™), with \m\ = mi + • ■ ■ + m n , and 
the rrij are the partitions of m characterizing a Young tableau. It then follows that 

(tiA) k = ]T Z m {A) (B.7) 

\m\=k 

in which the sum runs over all the Young tableaux with k boxes. We may now take for a 
generating function the hypergeometric function with matrix argument defined as 

p F q (a 1 ,...,a p ;P 1 ,...,P q ;t) = ±^ £ tf' ' 'th Z m (t). (B.8) 

fc=0 K - \m\=k vPl\m [Pq\m 



Then, one has 



F q (a 1 ,...,a p ;f3 1 ,...,(3 q ;d,b)= / p F q {a l ,...,a p ](5i,...,(5 qi gdg l b)dg 

•J (j 



IG 

oo 



^p J_ y^ [ Q i]m • • • K]m Z m (a)Z m (b) 

k=0 k\ | m | =fc [/3i] m • • • [f3 q ] m Z m (l) 



(B.9) 
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where 

n 

[«]r»=II(a-J + l)r«, (B.10) 

3=1 

(a) k = a(a + l)---(a + k-l) (B.ll) 

In the case of the unitary group G = U(n) the zonal polynomial Z m (A) is a Schur function, 
namely it is given by the Weyl formula for the characters of the representations of G 

i , / mj+n—j\ 

Z,M) = ^h^- (B.12) 

det(a i ) 

The Euler integral gives 

Zm{A) = Tn{ ® fk I Zm (rA)det(rr- n det(l - rf-^dr (B.13) 

r n (/3 - a)T n (a) [a\ m Jo<r<i 

where the integration is over Hermitian matrices r whose eigenvalues are between and 1. 
The Gamma function T n (a) is defined by 

n 

T n {a)=Tx n{n ~ 1 ^ 2 \{Y(a-i + l) (B.14) 

Using these notations, we have established the recurrence formula, 
iF q+ i(aci, ..., (Xp+i, /?i, ...,(3 q+ i] a, b) 



r n (/3 9 +i) 

p+i^q+iy^ii •••,«p+i, w, •••>/-vW'"> L V — : F7 vFTfl 7 



I det(r)^ +1 - n det(l - r)^ +1 ~ a " +1 - n p F q (a u .., a p ; p\, .., (3 q ; ra, b)dr 

J0<r<l 

(B.15) 



which proves (|B.5|) inductively. 



For the case of the confluent and Gaussian hypergeometric functions, \F\ and 2 -Pi, these 
formulae reduce to 

1 F 1 (a;/3;t) = - , ^ r / cire tr ^det(r) a -"det(l - rf — (B.16) 

1 n (ajl n {p — a) J0<r<l 

2 F 1 (a,(3;r,t) = - rw( ^ / rfrdet(r)/ 3 -det(l -r)^- 

r n (p)r„(7 - (3) Jo< r <i 

xdet{l -rty a (B.17) 

For the symplectic group Sp(n), we have similar formulae with hypergeometric functions. 
In the case p = 0, q = 0, and p = 1, q = 0, which is similar to (P-l| ) and ( |B.3| ), 
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D F (a,6) = / e tT{9a9 lb) dg 

JSp(n) 

det[2sinh(2aj6j)] 
A(o)A(6) 



:b.i8) 



1 F («; a, b) = f det(l - a^" 1 )-"^ 1 ^ 

= Ul - a ^ , A/L \ 1+a ' &J (B.19) 

A(o)A(6) 

where A(a) = n( a f ~~ a f) II&j- These formula are easily checked for the simplest case, n — 1, 
Sp(l), which is isomorphic to SU(2), with 

/ cos4>e i91 sm(pe ie2 \ . . 

P ~ I - sin 0e-^ cos 0e~^ J l j 

It may be useful to notice that these results for Sp(l) may also be derived immediately 
from the U(2) case by setting a — a,\ — —a 2 , and b — b± — —b 2 in (|B.1| ) and ( |B.3| ). Indeed 
the Lie group g of Sp(l) in ( B.20|) is derived from (p.4j ) with the condition 6* 3 = — 2 . The 



denominator A(a) for Sp(n) is also derived from the £/(n) case with the condition on the 
eigenvalues a = diag(ai, • ■ ■ ,a n , —ax, •• • , — a n ). This equivalence between U(n) and Sp{n) 
holds for n^l. Since the eigenvalues appear here in pairs (oj, — Oj), we have to replace 
n by 2n in the hypergeometric relations as shown in ( p. 19 ). For a general hypergeometric 
function, we have for the Spin) case, 



p F q (a 1 , ..., a p ; /3x, ..., (3 q ; a, b) 
/ p F q (a t , ..., a v \ fa, ...,P g ; agbg~ l )dg 

JSp{n) 

C 
— — — — — det[( p J r (? (ai - 2n + 1, ..., Op - 2n + 1; /9i - 2n + 1, 
A(a)A(o) 

/3 (? -2n + l;a i 6 i )) 2 

( p JF g (a; 1 — 2n + 1, ..., a p - 2n + 1; /?i — 2n + 1, ..., 
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P q - 2n + 1; -a^)) 2 ] 



(B.21) 
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